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ABSTRACT 

A recent high angular resolution extinction map toward the most opaque 
molecular globule, Globule 2, in the Coalsack Nebula revealed that it contains 
a strong central ring of dust column density. This ring represents a region of 
high density and pressure that is likely a transient and possibly turbulent struc- 
ture. Dynamical models suggest that the ring has formed as a result of a sudden 
increase in external pressure which is driving a compression wave into the Glob- 
ule. Here we combine the extinction measurements with a detailed study of the 
C^^O (1-0) molecular line profiles toward Globule 2 in order to investigate the 
overall kinematics and, in doing so, test this dynamical model. We find that the 
ring corresponds to an enhancement in the C^'^O non-thermal velocity disper- 
sion and non-thermal pressure. We observe a velocity gradient across the Globule 
that appears to trace two distinct systematic subsonic velocity flows that happen 
to converge within the ring. We suggest, therefore, that the ring has formed as 
two subsonic flows of turbulent gas merge within the Globule. The fact that the 
outer layers of the Globule appear stable against collapse yet there is no centrally 
condensed core, suggests that the Globule may be evolving from the outside in 
and has yet to stabilize, confirming its youth. 
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Introduction 



Understanding the initial conditions that lead to the formation of dense cores is impor- 
tant as it has direct implications for theories of star formation. Recent studies of starless 
cores reveal that the kinematics in the early stages of core collapse can be quite complex: in 
some cases the molecul ar tracers indicate simultaneously opposing motions (expansion and 
infall; iLada et al.ll2003l ). In the case of the starless Bok globule B68, the presence of both 
inward and outward motions, traced by molecular line e mission, has been interpreted to 



arise due to the oscillati on of the outer layers of the cloud (ILada et al.ll2003l : iRedman et al. 



20061 : iMaret et al.l 120071 ) . This idea is strengthened by modeling that suggests that the 



complex motions within starless cores can occur due to oscillating perturbations about a 
thermally-supported, pressure bounded structure (IKeto et al.ll2006l : iBroderick et al.l 120071 ). 
Indeed, recent numerical simulations reveal that small a mounts of rotation can cause initially 
unstable and collapsing cores to stabilize and oscillate (IMatsumoto fc Hanawall2003l ). 



The combination of high angular resolution extinction and molecular line observations 
toward cold, dense starless cores can trace the properties and kinematics of the dense gas and 
allow us to gain insight into the physical nature of core and, thus, star formation. Extinction 
measurements are particularly useful for tracing the internal structure of cores because they 
probe the material along the line-of-sight and are not susceptible to optical depth effects (as 
in the case of ^^CO molecular line emission) or assumptions about the dust temperature or 
emissivity (as in the case of millimeter continuum emission). Thus, the combination of dust 
extinction observations with optically thin molecular line observations can accurately reveal 
the masses and densities of the dust as well as the kinematics and dynamics of the gas. 

In this paper we combine extinction observations with a detailed kinematic study of the 
molecular line profiles observed toward Tapia's Globule 2 in the Coalsack Nebula. Because 
Globule 2 does not contain a central protostar, is one of the densest cores within the Coalsack 
complex, and is thought to be in a very early stage of condensing, a detailed study of 
its kinematic structure could provide important details on the processes leading to core 
formation. 

The Coalsack Nebula is a nearby dark molecular cloud complex that extends almost 6° 
on the sky (~ 15 pc at its distance of 150 pc: ICambresvlll999) and appears complete ly devoid 



of star formation activity (jNyman et al. 1989 : Bourke et al 



19951 : iKato et al.lll999[ ). It has a 
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total mass of ~ 3550 Mq and is characterized by complex, filamentary molecu lar structures 



which contain many dark cores with kinetic temperatures, T^^, of 8-11 K (INvman et al. 



1989). Th e darkest and densest of these cores are Tapia's Globules 1,2, and 3 (ITapia 
Boklll977h . 



1973 



Using deep near-infrared imaging observations obtaine d with the Europ ean Southern 
Observatory (ESO) 3.6m New Technology Telescope (NTT). lLada et al.l (120041 ) were able to 
examine the structu re of Globule 2 in significantly grea ter detail than previously possible. 
The NICE method (ILada et al.l Il994j : lAlves et al.l Il998l ) was used to generate a visual ex- 
tinction (Ay) map that was constructed by measuring the individual infrared extinctions to 
approximately 24,000 stars located behind Globule 2. This resulted in an extinction map 
across the cor e with significantly improved angular resolu tion (~ 15") compared to previous 
maps (~ 60": ljones etaDll98ol ^983 : iRacca et al.ll2002h . 



The extinction map revealed Globule 2 to have a mass of ~6.1 Mq, a mean density of 
■^3xl0^cm~^ and Jeans mass, Mj, of ~10.4 Mq. From near-infrared polarimetry obser- 



vations, IJones et al.l (119841 ) suggested that Globule 2 is supported either by turbulence or a 
magnetic field and appears stable. This is also consistent with the fact that the derived mass 
is less than Mj suggesting it is unlikely to fragment further. Moreover, a Bonnor-Ebert (BE) 
fit to the outer regions of the Globule's radial extinction profile suggests that it is globally 
stable against collapse (ILada et al.ll2004l ). The Globule's derived BE mass is ~ 6 Mq. Even 
though this is comparable to the derived mass calculated from the dust, the high angular 
resolution infrared extinction map toward the Globule does not reveal a centrally condensed 
core as would be expected if it was a true BE sphere. Instead, the most prominent feature 
within the extinction map is a strong central ring of dust column density. This ring likely 
represents a trans ient structure that has formed from a region of high density and pressure 
( iLada et al.ll2004l ) . Because there is no source of internal pressure to support the ring against 
gravitational collapse, it is thought that the ring is not in dynamical equilibrium with its 
surroundings and, given the short timescales for inward collapse (2x10^ years), most likely 
represents an extremely early stage of evolution. 



From a single C^^O (2-1) spectrum toward the center of the globule, iLada et al.l (120041 ) 
noted a double-peaked profile comprising a bright narrow line and a fainter broad line. 
Based on the abse r ice of this double-lined profile in low er angular resolution spectra (e.g.. 
Brooks et al.l 1 19761 : iNyman et al.l Il989l : iKato et al.l Il999l ) , it was suggested that the narrow 
line component most likely arises within a small region an d that the broad line component 
traces a more spatially extended feature (ILada et al.ll2004l ). The hne width of each blended 
component was found to be super-thermal indicating that significant non-thermal, or turbu- 
lent, motions most likely dominate. Dynamical models for Globule 2, however, suggest that 
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the ring structure has formed as a result of a sudden increase in ext ernal pressure which is 
driving a compression wave into the Globule (IHennebelle et al.ll2006l ). 



Using the Mopra 22-m telescope, we have obtained the first complete high angular and 
spectral resolution molecular line map of Globule 2. Because ^^CO is optically thick toward 
most dense cores, and is thus only tracing the outer layers, we choose to conduct this study 
using the optically thin C^^O (1-0) transition which is far superior in tracing the dense gas 
associated with the central core. In combination with the extinction map, these data should 
allow us to test the various formation mechanisms for the ring feature; is it a h igh-density 
transient structure formed via the dissipation of turbulent gas (ILada et al.l 120041) or does it 
trace the propagation of a compression wave through the Globule (IHennebelle et al.ll2006l )? 
Moreover, these data will enable the relative spatial distributions of the two molecular line 
components to be definitively traced. In particular, does the narrow-line component originate 
exclusively in the dense ring? Does the broad-line component primarily trace the turbulent 
envelope of the Globule? Answers to these questions will directly test the idea that Globule 
2 is in perhaps the earliest stage of evolution: a dense core in the process of formation. 



2. Observations 

Molecular line maps of Globule 2 in the Coalsack Nebula were obtained with the Mopra 
22-m telescope in 2006 July. The 8 GHz spectrometer MOPS was used to simultaneously 
observe a number of transitions including C^'^0 (1-0)0. The spectrometer was used in 'zoom' 
mode such that one spectral window covered each of the lines of interest. Each window 
was 137.5 MHz wide and contained 4096 channels in both orthogonal polarizations. This 
produced a velocity resolut ion of 0.09 kms~ ^ for the C^^O transition. At these frequencies 



the Mopra beam is ~ 33" flLadd et al.ll2005f ) 



The observations were obtained using the On-The-Fly (OTF) mode by scanning in 
both Right Ascension and Declination. All maps were Nyquist sampled. A fixed off-source 
position was used for all maps (a=12:29:57.1, 5=— 63:32:21.8, J2000). The nearby SiO maser, 
R Car, was used to check the pointing accuracy approximately every hour and was typically 
better than 9". The system temperatures were checked frequently using a paddle and were 
typically 180 K. 

A total of fourteen 5'x5' OTF maps were combined (using the system temperature 



^The other transitions observed included ^^CO, ^^CO, HCO+, CH3OH, and CS and will be presented in 
future papers (Saul et al., Butner et al. in prep.). 
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weighted median) to produce the final map. The data were initially baseline subtracted and 
converted to a temperature scale using the AIPS^^ task livedata. All raw data are in the 
scale. To c onver t to the main beam brightness, T^^, we use the beam efficiencies listed 



T 



4, 
in 



Ladd et all fl2005[ ). 



The OTP maps were combined and gridded into data cubes using the AIPS"'"^ task 
gridzilla. To produce the final map, a Gaussian smoothing kernel with a PWHM of 0.55' 
was applied to the data. This resulted in the final cube having an effective angular resolution 
of ~ 46". The final C^^O map has a T^ rms noise sensitivity of 0.08 K channel"^. 



3. Results 

Figure [1] shows the visual extinction map of Globule 2 (15" angular resolution). This 
map was generated from de ep near-infrared images by measuring the extinction toward 



~ 24,000 background stars (ILada et al.l l2004j ). The ring is clearly seen as an extinction 
enhancement surrounding a central depression. Overlaid on this map is the C^^O (1-0) 
integrated emission obtained with the Mopra 22-m telescope (the emission was integrated 
over the velocity range of —7.0 to —5.0 kms~^; 46" angular resolution). These data show 
that the C^^O integrated emission follows the overall morphology of the extinction. There 
are noticeable differences, however. The C^^O integrated emission does not trace completely 
the extinction ring. In particular, it appears that the extinction peaks toward the eastern 
arc of the ring, whereas the C^^O integrated emission peaks toward western arc. Some of 
these differences may arise due to the different angular resolutions of the two datasets (15" 
and 45" respectively). When smoothed to an angular resolution of 45", the extinction image 
shows a much shallower contrast between the ring and the central hole. 

Pigure [2] shows the C^^O (1-0) channel maps across the Globule overlaid with contours 
of the smoothed visual extinction image (the Ay map of Pig. [1] was smoothed to 46" to 
match the angular resolution of the C^^O data). A clear velocity gradient is evident across 
the Globule. 



3.1. Gaussian component analysis 

The individual C^^O (1-0) spectra across the Globule reveal that the emission pro- 
file differs quite considerably on and off the extinction ring. Based on the detection of a 
'narrow-line' (v2,s'ij~— 5.64 kms~^) and 'broad-line' (v2,5i?~— 5.83 kms~^) component within 
the C^^O (2-1) spectrum toward the ring and the absence of the narrow-line component 
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within lower-angular resolution surveys, it was suggested that the broad line was spatially 
extended across t he Globule while the narrow line was restricted to the extinction ring 



( iLada et al.ll2004j ). Thus, to characterize and separate these two emission features and de- 
termine any spatial differences between them, we fit Gaussian profiles to each spectrum 
within the map. 

Because we expect the broad-line component to contribute to the observed spectra at all 
points across the map, but the narrow-line component to only those spectra coincident with 
the extinction ring, we attempt to fit both a one- and two-component Gaussian profile to 
each spectrum within the map. Initial guesses for the vlsr, hne widt h, and peak t emper ature 



for the Gaussian profiles were based on the observed profiles from iLada et al.l (120041 ) . All 
were free parameter^. For example, when fitting a one-component Gaussian profile we use 
as input the values determined for the 'broad' line. When fitting a two-component Gaussian 
profile we input the values derived for both the 'broad' and 'narrow' lines. At each position 
within the map we select either the one- or two-component Gaussian fit based on which had 
the lower residuals. We then assigned each fitted Gaussian profile to either the 'narrow' or 
'broad' line component based on its derived vlsr- We find that ~ 70% of the spectra were 
well fit by a one-component Gaussian profile, while the remaining ~ 30% were best fit by 
two components. 

As expected, we find that the spectra toward the extinction ring are best matched by a 
two-component Gaussian profile, while the spectra external to the ring are best matched by 
a single component. Contrary to the original prediction, however, the central vlsr derived 
for the emission external to the ring does not correspond to a single vlsr, nor does it simply 
correspond to the 'broad' component. Instead, the vlsr shows a range of values (from 
—6.0 to — 5.7kms~^) which suggests that both the 'broad' and 'narrow' lines are located 
external to the ring. It appears, therefore, that the narrow-line component may not simply 
arise within the ring, but instead may correspond to a spatially separate velocity component 
extended across the Globule. 

This can be seen in the C^^O channel maps of Figure[2) the component at ~— 5.64 kms~^ 
peaks in the ring, but also extends beyond it. The broad-line component (~— 5.82 kms^^) 
also arises both in and out of the ring, the only difference being it appears to trace the 
molecular material on the opposing side of the extinction ring to the narrow-line compo- 
nent. Figure [3] shows the measured peak temperature (Tmb), central velocity (yLSR), and 
one-dimensional velocity dispersion (a) for the two components from the Gaussian fitting 
procedure. We find that the individual measured central vlsr and a for the two components 



^the Gaussian fitting was performed within IDL using the gaussfits procedure. 
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show a dispersion in their values. This suggests that the kinematics of the region are quite 
complex and not well matched simply by two Gaussian vlsr components toward the ring 
and a single vlsr component external to it. 

The absence of the narrow-line component in previous molecular line surveys is not 
surprising given their limited angular and velocity resolutions. Taking an average C^^O (1- 
0) spectra oyer t he region contained within the beam of the H2CO observations (4'4; 



Brooks et al.lll976l ). we find that the molecular line profile is characterized by a single emis- 



sion line, with vlsr^ —5.8 kms ^ and line width, AV, of 0.6 kms ^. 



3.2. Moment analysis 

Rather than refer to the individual components as the 'narrow' and 'broad' lines based 
on their original vlsr. and Gaussian identification, these data reveal that the C^^O emission 
displays complex velocity features across the Globule. Thus, in this section we characterize 
the emission via a moment analysis which is independent of Gaussian deconvolution and 
assumptions about individual components. 

Over the extent of the Globule, the combination of these velocity features manifest 
themselves as a velocity gradient. This velocity gradient is seen in the C^^O channel maps 
(Fig. [2]). While the brightest emission coincides with the extinction ring, we also see the red- 
and blue-shifted emission tending to arise on opposite sides of the extinction ring. This can 
also be seen in the position- velocity, (i-V) , diagrams of Figure HI The left panel shows the 
(i-V) diagram averaged over the extinction ring in Declination (the emission was averaged 
over A6 of 6' centered on the ring). The right panel shows the Right Ascension-averaged 
diagram (the emission was averaged over Aa of 18' centered on the ring). The velocity 
gradient is more apparent in the Declination-averaged {i-V) diagram. 

The zeroth moment map, or integrated intensity image, is defined as 




where Tmbi^) is the main beam brightness at a given velocity v. The image was generated 
by integrating all of the emission over the range of —7.0 to —5.0 kms~^ and is shown in 
Figure [5] overlaid with contours of the smoothed visual extinction image. The extinction 
ring is clearly seen as an enhancement in C^^O emission. 

The first moment map. Mi, is a measure of the intensity (in this case Tmb) weighted 
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velocity field and is defined as 



Ml 



/ Tmb{v)vdv 
J Tmbiv) dv 



(ISault et al.lll995l ). The C^^O first moment map is also shown in Figure [5] and clearly shows 
the aforementioned velocity gradient. Toward the edges of the extinction ring the velocity 
field ranges from —6.0 to —5.7 kms~^ Toward the peak of the extinction ring, the vlsr of 
the emission is ~ —5.8 kms^^. 

The intensity weighted velocity dispersion is represented via the second moment, M2, 
which is defined as 



l jTmbiv)iv-M'i 



'■ dv 



f Tmb{v) dv 



The C^^O second moment map is also shown in Figure [51 We see a clear enhancement in 
the velocity dispersion toward the extinction ring. Moreover, the morphology of the velocity 
dispersion and the Ay appear to match very well. Indeed, of all the measured quantities, 
the velocity dispersion matches best with the visual extinction. 

In addition to the velocity gradient, the line profiles clearly show an asymmetry which 
changes across the Globule. To quantify the asymmetry of the emission profiles we have 
calculated the 'skewness' of the emission at each position within the map. We follow 



Gregersen et al.l (119971 ) and define skewness as a dimensionless ratio of the third moment 
to the 3/2 power of the second moment where both are normalized by the first moment, 
such that 



, / Tmb{v){v - vlsr)'"^ dv 

skewness = / 

j Tmb{v) dv 



f Tmb{v){v - VLSRfdv 

J Tmb{v) dv 



3/2 



where vlsr is the mean velocity of the Globule (~— 5.84 kms~^). Figure [5] also shows the 
skewness image and reveals that, toward the west of the extinction ring, the profiles show 
a negative (blue) asymmetry, while toward the east, the profiles typically show a positive 
(red) asymmetry. 



3.3. Temperature and density tracers 



In addition to the C^^^O map, we have also obtained ^^CO and CS maps toward the 
Globule. While these results will be discussed in detail in Saul et al. (in prep), we include 
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them here to help with the interpretation of the C^^O resuhs. Because ^^CO is typically 
optically thick, it can be used to trace the kinetic temperature of the gas. Because it has 
a large dipole moment, CS has a high critical density, and is typically excited coUisionally 
only in high-density gas (n>10^ cm~^). Thus, the combination of ^^CO and CS can trace 
the warm, dense gas within the Globule. 

Figure [6] shows the kinetic temperature map overlaid with contours of the smoothed 
visual extinction map. We find that the kinetic temperature displays a range of values 
over the Globule (from ~ 8 to lOK). The highest temperatures are not associated with the 
extinction ring, but are located to the north-west (Tk of 9.6 ± 0.5 K compared to ~ 7.6 ± 
0.4 K in the south-east). The kinetic temperature of the gas associated with the extinction 
ring is ~ 8.8 ± 0.4 K. 

We find that the dense gas (as traced by the CS integrated intensity) is associated with 
the extinction ring (Fig. [7]). While the morphology of the CS emission and extinction do not 
match exactly, there is a general correspondence between the two. However, the morphology 
of the CS emission agrees very well with the C^^O emission: both peak on the western arc 
of the extinction ring. The fact that the C^^O and CS emission match so well suggests that, 
chemically, this core does not appear to be evolved. Because carbon-chain molecules are 
considered 'early-time' molecules, they are abundant in chemically young cores and depleted 
in more evolved, star-forming cores. This is in contrast to molecules such as N2H+ and 
NH3, which ar e considered 'late-time' molecules and, thus, are more abundan t in chemically 
evolved cores (IBergin fc Langerlll997l : lAikawa et al.ll2005l : ISuzuki et al.lll992l ). 



4. Discussion 

4.1. Comparison of the dust extinction to molecular abundances 

We find that the morphology of the C^^O (1-0) integrated intensity generally matches 
the visual extinction. By smoothing the Ay image to match the angular resolution of the 
q18q jj^a^p^ f^gj^ make a point-by-point comparison between the integrated C^^O emission 
and the visual extinction across the Globule. 

Figure [8] shows the C^^O integrated intensity as a function of Ay- The individual 
data points are shown as small dots while the filled circles trace the overall trend and were 
generated by taking the median of the integrated intensities within evenly-spaced Ay bins. 
Each Ay bin is 0.7 mag wide; the errors bars in Ay represent the range in each bin. The 
error bars in the integrated intensity represent the dispersion in the data. We find that the 
C^^O integrated intensity is generally correlated with visual extinction across the complete 
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range of measured extinction values. However, for any given Ay, there appears to be a 
range in the measured integrated intensities. Similar to what is observed toward other dark 
clouds, we find that C^^O appears to be depleted only at the highest measured extinction 
(Ay ^ 9.5 mag). Figure [S] also shows a similar plot for the CS integrated intensity as a 
function of Ay and reveals that CS is also only depleted at the very highest extinctions. A 
simil ar analysis of the integrated intensity as a function of Ay toward the denser core 



B68 flBergin et al.ll2002l ) shows that C^^O is depleted for Ay > 10 mags Thus, we do not 
see significant molecular depletion toward the Globule primarily because its mean column 
density is too low. The lack of a high density, centrally condensed core further suggests that 
the Globule is relatively young. 

A least-squares fit to the C^^O data (solid line) for 4 mag<Ay <9.5 mag reveals that they 
are well matched to a linear fit with a constant slope over the complete range of measured 
extinction values. The derived linear fit is 

/(Ci^O) = (0.17 ± 0.01)A - 0.04 ± 0.07 K kms'^ 

The relative abundance of C^^O to molecular hydrogen (H2) can be estimated from 
comparisons between the derived values of C^^^O column density, N(C^^O), and Ay. A 
least-squares linear fit to these data also show a linear correlation which is given by 



N{C^^O) = (1.4 ± 0.1) X 10^M„ + (0.8 ± 0.4) x 10 



a4 -2 
^ cm 



Our derived slope is slightly lower than that calculated by iLada et al.l (Il994l ) within IC 
5146 (~2.3 xlO^^ cm~^) but is c onsistent with the range observed toward oth er dark clouds 
(0.7-3.5 X 10^^ cm~^ mag~^; see lHayakawa et al.lll999l : lHarjunpaa et al.ll2004l and references 
therein). The tight correlation in the integrated intensity and, hence, C^^O column density 
over all measured Ay, suggests that there is little variation in the molecular abundance of 
C^^O as a function of Ay and that the assumption of a single excitation temperature across 
the Globule is most likely valid. 

By converting the measured Ay to a H2 column density, where N(H2)/Ay=10^^ cm~^ mag" 
we can also derive the relation between the C^^O and H2 column densities. A least-squares 
fit to plots of the data measured as N(C^^O) and N(H2) gives the relation 



cm-2 



N{C^^O) = (1.4 ± 0.6) X W'^N{H2) - (0.8 ± 0.4) x 10 
'he d erived slope is comparable to the value of ~ 1.8 x 10~^ calculated by iFrerking et al. 



(119821 ) toward the dense cores within Taur us and p Ophiuchus and to the valu es determined 



for the dense cores L 977 and FeSt 1-457 flAlves et al.lll999l : lAguti et all 120071 ) 
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Even though we see very httle overall C^^O depletion toward the Globule, the apparent 
differences between the peaks of the visual extinction and C^^O in the ring may arise as 
a result of depletion. Because the western arc of the ring has a higher mean extinction 
compared to the eastern arc, it is possible that the C^^O is more depleted toward the western 
arc. These difference are seen in Figure [9] which plots the mean extinction, C^^O integrated 
intensity, and the ratio of C^^O integrated intensity to Ay across the Globule. For this 
analysis, the Ay image was smoothed to match the angular resolution of the C^^O. These 
plots were generated by taking the mean values in a 45" strip of constant Declination across 
the Globule. The enhanced extinction toward the ring is clearly seen, along with its central 
hole which is evident as an ~ 10% decrement in the column density (top panel). This is in 
contrast to the C^^O integrated intensity (middle panel) which does not show a significant 
decrement. These differences are most obvious in the ratio of the integrated intensity 

to Ay (lower panel). We find that the ratio is essentially constant toward the western arc 
of the ring (positive offsets), but drops by ~ 20% toward the eastern arc (negative offsets). 
This is consistent with the fa ct that depletion is typically only seen toward regions with the 
highest Ay (> 10 mags; e.g.. lBergin et al.ll2002l ). Because the mean Ay within most of the 
extinction ring is very close to this threshold, it may be that we are seeing the early onset of 
detectable C^^O depletion in the western arc of the extinction ring. It is possible that more 
depletion is seen here simply because it has a higher column density. 

This is consistent with a recer it study of the solid- state H2O ice column densities toward 
the Globules within the Coalsack (IRodgers et al.ll2007l ). Toward the outer edges of Globule 2 
they find only a small amount of oxygen depletion. Because the time required to accrete the 
observed amount of oxygen is significantly shorter than the tim e for full oxygen dep letion, 
they conclude that the Globule is relatively young (< 10^ years; iRodgers et al.l 120071 ). The 
fact that this dep letion time scale is comparable to the Globule's dynamical time scale 
( iLada et al.l l2004j ) and the small amount of detectable depletion confirms the youth of the 
Globule. 



4.2. Global properties of the Globule 



While Globule 2 is one of the densest cores within the Coalsack Nebula, with a mean den- 
sity nfH^) of ~2.7xl 0^ cm~^, it is less dense than typical cores found in Taurus (~ 1.2 x lO'^cm"^; 



Ladd fc CovevI l2000h and the Pipe Nebula (^7.3x10^ cm'^; Ihada et al.l 120081 ). Because it 



is quiescent and does not harbor a central protostar, it is considered a dense, starless core. 
Moreover, the high Jeans mass to total mass ratio and the Bonnor-Ebert fit to the density 
profile for the outer regions of the Globule, suggest that the Globule is stable against col- 
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lapse (ILada et al.ll2004l ). However, the discovery of an extinction ring within the Globule 
has lead to speculations that it may be a dynamically unstable, transient core on the verge 
of condensing to form a protostar. With our C^^^O map, we can now test this idea and 
determine whether the Globule is gravitationally bound, in virial equilibrium, magnetically 
super-critical, rotating, thermally dominated, and/or pressure confined. These properties 
can reveal important clues to the nature and evolutionary status of the Globule. 



4.2.1. Gravitationally bound? 



From the single line profile toward Globule 2, iLada et al.l (l2004j ) suggested that the 
Globule is gravitationally bound. We now extend this analysis to include all the data mea- 
sured toward the Globule. By comparing the escape velocity, Vesc, to the average measured 
three-dimensional velocity dispersion across the Globule, we can determine if the Globule is 
gravitationally bound. The escape velocity is given by 

Vesc = V^GM/R 

where M and R are the Globule's mass and radius respectively. We use the values estimated 
from the dust extinction for the mass (M=6.1 M(7-, ) and the radius calculated from the 



Bonnor-Ebert fit (R=290", 0.21 pc, or 6.51 xlO^^ cm; iLada et al.ll2004f ). Given these values, 
we estimate that the Vesc for Globule 2 is 0.50 kms~^. 

The average three-dimensional velocity dispersion, a^d, is calculated via the expression 

oSd = V 3a^ + Salt 



where a is the one-dimensional isothermal sound speed (0.19 kms"^ for hydrogen in a 10 K 
gas) and ant is the one-dimensional non-thermal velocity dispersion, i.e. ant = a/o"^ — afj^. 
The measured velocity dispersion, a, was estimated from the average in the second mo- 
ment map across the Globule and is ~ 0.18 kms~^. Assuming the gas is at 10 K, the 
one-dimensional thermal velocity dispersion, ath, is 0.05 kms^^. Thus, the non-thermal 
velocity dispersion, ant, is ~ 0.17 kms^^ and implies that the turbulent motions are sub- 
sonic. Using these values we find that the average three-dimensional velocity dispersion, a^d, 
is 0.44 kms~^. Because a^d is formally less than the escape velocity, it is likely that the 
Globule is indeed gravitationally bound, although only marginally so. 



The one-dimensional virial velocity dispersion, amriai, was also calculated via the ex- 
pression 
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IGM 

and was found to be 0.16 kms~^. This is lower than both the sound speed in a 10 K gas 
(0.19 kms~^) and the measured C^^O velocity dispersion across the Globule (0.18 kms^^). 
The fact that the dmriai is lower than the measured velocity dispersions indicates that the 
Globule is not virialized. Thus, it appears that the Globule is only marginally bound and 
not in virial equilibrium, implying that is it extremely young. 



4-2.2. Magnetically super- critical? 



The role of magnetic support within the Globule can be measured via the ratio of the 
observed mass to magnetic flux ($) represented in terms of the critical value. We follow 
Crutcherl (120071 ) and define this ratio as 



A 



(M/*) 



observed 



(M/$) 



7.6 X 10 



-21 



critical 



N{H2 
B 



where N(H2) is in cm ^ and B is in yuG. To calculate the magnetic field, B, we use the 
Chandrasekhar &: Fermil (119531 ) equation 



B 



O-0 



10 



2.7 X 103 



jj,G 



where is the one-dimensional velocity dispersion and c tb is the dispersion in the measured 
polarization angles. Using (T^=0.18 kms~^, (Te=0.66 rad (I Jones et al.lll984l ). and a molecular 
hydrogen density, n(H2), of 2.7x 10^ cm"^, we find that B~10 /xG. Thus, the directly observe d 
ratio, Xobserved, IS ~ 4.8. Applying corrections for geometrical effects of 1/3 (ICrutcherl 120071 ). 
we find that A is ~ 1.4. A value of A greater than unity implies that the magnetic field is 
super-critical. Thus, we find that the Globule is marginally super-critical which suggests 
that the Globule is not dominated by magnetic support and that it is likely bound if not 
already in the process of gravitational contraction. 



4.2.3. Rotating? 



As mentioned previously, the channel maps, {(!--V) diagram, and first moment map (the 
intensity weighted velocity field) all show signs of a velocity gradient across the Globule. In 
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the C^^O first moment map (Fig. [5]), the lines of constant velocity are nearly parallel and 
are aligned north-south. If the velocity gradient is linear when projected on the sky, we 
can relate the observed velocity of a ny point on the projected surface of the Globule to the 
velocity gradient via the expression (IGoodman et al.lll993l : iLada et al.ll2003l ) 



vlsr 



dv . „ dv . ^ . ^ 
f + — Aa cost^ + — Ao smt^ 
ds ds 



where Aa and A5 are the offsets in Right Ascension and Declination (in arcsec), dv/ds is the 
magnitude of the velocity gradient in the plane of the sky, 9 is its direction (measured from 
east of north), and vq is the systemic radial velocity of the Globule. The magnitude and di- 
rection of the velocity gradient and the Globule's systemic velocity were determined through 
a least-squares fit to the observed C^^O velocity distribution (as shown in Fig. [5]). We find 
that dv/ds = 1.00 ± 0.01 kms-^pc-^ 6 = 90.1 ± 0.3°, and = -5.81 ± 0.01 kms'^ This 
vel ocity gradient is v ery similar to that found in the C^^O survey of the Coalsack Nebula 
by lKato et al.l (119991 ). Even though these data were under-sampled {2'.7 beam on a 4' grid), 
they reveal a velocity gradient across Globule 2 of ~ 0.90 kms^^pc^^. If the Globule was 
undergoing solid-body rotation whose axis is perpendicular to the line-of-sigh t, then this 
veloc ity gradient would correspond to an ang ular velocity of 2.9 x 10"^^ s"^ flKato et al.l 
1999h . 



It is possible that the velocity gradient arises due to the rotation of the Globule. As- 
suming solid-body rotation, the ratio of rotational kinetic energy to gravitational potential 
energy, (3, can be used to determine th e signi ficance of rotation to the overall djTiamics of 
the Globule. Following iGoodman et al.l (Il993l ). we calculate (3 via the expression 



P 



(l/2)/w;2 Ipw^R^ 



qGMyR 2q GM 



where M and R are the mass and radius of the Globule respectively, the moment of inertia 
is given by / = pMR"^, q is defined such that the gravitational potential energy is qGM'^/R, 
and w = {dv / ds)/ sin i. We assume a sphere with a r~^ density profile, which implies p/q is 
0.22, and that sini=l. Using M=6.1 Mq, R=6.51x10^^ cm, and dv/ds=i m kms~^pc~ 



we firi d that /3=0.04. This value of (3 is similar to other cores: for example, IGoodman et al. 
(I1993I ) find that most clouds have /? < 0.05. 



Thus, we calculate that the apparent rotation does not contribute significantly to the 
overall stability of Globule 2. In any event, it appears that the observed dynamics of the 
Globule are more complicated than simply solid body rotation and it is likely, therefore, that 
the observed velocity gradient is not tracing a coherent rotating structure. Nevertheless, the 
energy in the velocity gradient is unimportant relative to gravity. 
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4- 2.4- Thermally dominated? 

The one- dimensional C^'^O velocity dispersions {a ~ 0.18 kms^^) are ~ 3.6 times 
broader than what would be expected if the profiles were purely dominated by thermal 
motions. This suggests that significant non-thermal motions characterize the velocity field 
of the gas and that turbulence may play a role in the internal pressure of the Globule. To 
quantify the degree to which turbulence may be contributing to the internal pressure, we 
calculate the ratio of the thermal to non-thermal (turbulent) pressure as given by 



R„ 



a, 



nt 



where a is the one-dimensional isothermal sound speed (0.19 kms~^ for hydrogen in a 10 K 
gas) and a^it is the one-dimensional non-thermal velocity dispersion. Taking the average 
ant across the Globule, we estimate the ratio of thermal to non-thermal pressure, Rp, is 
~ 1.2. Thus, it appears that thermal pressure is comparable to, or slightly greater than, the 
non-thermal, or turbulent, pressure within the Globule. 

While it is useful to calculate average values across the Globule, there are clear dif- 
ferences in the C^^O emission within and either side of the extinction ring. Indeed, if we 
take the typical measured velocity dispersion across the extinction ring and compare it to 
two regions on opposing sides of the ring, we find that the velocity dispersion toward the 
ring (cr of 0.20 ± 0.01 kms~^) is ~ 1.3 times greater than that observed outside the ring 
(0.15± 0.01 and 0.17± 0.01 kms~^ to the east and west of the ring respectively). Using 
these values we find that the non-thermal velocity dispersions (cr„t) in and out of the ring 
are 0.19 and ~0.15 kms~^ respectively. Thus, we see turbulent sonic motions within the 
ring and turbulent subsonic motions external to it. 

The differences between these regions are most obvious when we consider the ratio of the 
thermal to turbulent pressures, Rp. within the ring, i?p is ~ 1, while external to the ring Rp 
is ~1.6. This suggests that thermal pressure is more dominant external to the ring, whereas 
toward the ring, the contributions from thermal and turbulent pressure are comparable. 

This can also be seen in the left panels of Figure [TO] which shows the C^^O non-thermal 
velocity dispersion {ant) and the ratio of thermal to non-thermal pressure (Rp) as a function 
of Ay for all positions in the map. In the upper panel the dotted horizontal line marks the 
one-dimensional isothermal sound speed. In the lower panel the dotted horizontal line marks 
the point at which thermal and turbulent pressure support are equal. The dashed vertical 
line marks the approximate visual extinction at the edge of the extinction ring (~ 9.5 mags). 
These plots show that the general trend is for higher non-thermal velocity dispersions, cr„t. 
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and lower values for Rp as Ay increases. Coincident with the extinction ring (Jnt transitions 
from subsonic to trans-sonic and Rp transitions from greater than to less than 1. Thus, it 
appears that the outer regions of the Globule are thermally dominated, while toward the 
extinction ring non-thermal pressure may be more important. 



4-2.5. Pressure confined^ 



To further characterize its nature, we have also calculated the internal gas pressures 
within the Globule. The total pressure within the Globule is measured by a combination 
of the pressures due to thermal (Pth) and non-thermal motions (P„t). The thermal and 
non-thermal pressures were calculated via the expressions 



Pi 



th 



pc? and Pnt = po".^ 



nt 



where p is the mean density, p = 3M/{AttR^), a is the one-dimensional isothermal sound 
speed, and ant is the one-dimensional non-thermal velocity dispersion. We find that Pth/k ~ 
2.8 X 10^ K cm~^ and Pnt/k ~ 2.2 x 10^ K cm~^. Thus, the total pressure, Ptotai/k, due to 
thermal and non-thermal motions is ~ 5.0x10^ Kcm~^. 

Because the Globule lies within the larger Coalsack molecular cloud complex, it is 
possible that the weight of the Coalsack molecular cloud may also be a significant source of 
confining pressure. The pressure due to the weight of the cloud, Pdoud, is given by 



P. 



cloud 



3tt 
20' 



5G 



where G is the Gravitational constant, S is the mean mass surface density of the cloud 
(E = M/vri?^), and (t> n is a dimensionless corr ection factor to account for the non-spherical 
geometry of a cloud (IBertoldi fc McKeelll992r). Given the Coalsack molecular cloud has a 
mass of ~3550 Mq, radius of ~ 7.1 pc (INyman et al.lll989l ). and assuming 0g = 1-6, we find 
that Pcioud = 0.8 X 10^ K cm'^. 

This pressure is more than a factor of 5 lower than the internal pressure calculated from 
a combination of the thermal and non-thermal motions which suggests that it is unlikely that 
the weight of the Coalsack molecular cloud contributes significantly to the overall confining 
pressure of the Globule. Interestingly, this pressure is comparable to the pressu re resulting 
from the interstellar medium, Pisivi/k ~ lO'^Kcm"^ ( IBertoldi fc McKedll992l ). Thus, we 
find that the total internal pressure is greater than both PisAi/k and Pdoud/ k and that the 
Globule is super-critical. This adds to earlier suggestions that the Globule is gravitationally 
bound. This is in contrast to dense cores found within the Pipe Nebula most of which are 
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thought to be the rmally dominated, unbound cores that are confined by external pressure 
( iLada et al.ll2008l ) and suggests that the Globule is bound by gravity. 



4.3. The Origin of the Ring 

4-3.1. External pressure driving a compression wave? 



Recent dynamical models for Globule 2 (IHennebelle et al.ll2006l ) suggest that the ring 
structure has formed due to a rapid increase in the external pressure on a core with very 
little initial turbulence. In this picture, the ring structure traces the dense regions swept up 
by the compression wave as it propagates through the core. The central depression exists 
be cause the compres s ion Av ave has not yet converged in the center. The models generated 
by IHennebelle et al.l (120061 ) predict that toward the center of the Globule the line profiles 
should be double-peaked, while toward the outer edges of the Globule the profiles should 
show a single line. 

With our C^^O m ap we can now te s t this model. Due to the receding and approaching 
compression wave, the IHennebelle et al.l (120061 ) model predicts that positions within 50" of 
the center should show a double-peaked line profil e. Figure [11] shows a comparison between 
our data and the model of IHennebelle et al.l (120061 ) for different impact parameters, b, across 
the Globule. The spectra were generated by averaging over annuli at the impact parameters 
6=0, 0.04, 0.08 and 0.1 2 pc. The example model spectra were estimated from figure 2 of 
Hennebelle et al.l ( l2006l ) and are plotted so that the center velocity (their V=0.0 kms~^) 
corresponds to the v^^^^ of the Globule (—5.8 kms~^). The individual line centers, widths 
and relative intensities were taken from their figure, while the peak temperatures were scaled 
to approximately match the data. 

From this comparison we find that at larger impact parameters (6=0.08 and 0.12 pc) 
the models match the data well. These regions correspond to the edges of the Globule, 
off the extinction ring. Toward the ring, however, we find the models differ considerably 
from the data. While we do observe a double-line profile at impact parameters b < 0.04 pc, 
the separation of the line centers and the relative intensities of the two components are no- 
ticeably different from what the models predict. Toward the center of the ring (6=0.0 pc), 
we find that the two components are separated by ~ 0.24 kms~^, with the red component 
dominating. Here the models predict the line profiles should be separated by 1 kms~^ and 
have the same relative intensity. Similarly, at an impact parameter 6=0.04 pc (which com- 
pletely encompasses the extinction ring), we also see a dominant red component separated 
by ~ 0.21 kms~^ from the fainter blue component. 
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The discrepancies between the observed hne profiles and the model predictions could 
arise simply because of the assumptions of a symmetric compression wave. The apparent 
fragmentation within, and the non-symmetrical nature of, the extinction ring suggests that 
the mechanism that has given rise to the ring is unlikely to be smooth and symmetric. 
Moreover, the predicted line separations could be higher than that measured due to an 
overestimation of the speed at which the material is imploding within the models. While the 
model predicts the integrated intensity and velocity dispersion should increase within the 
extinction ring, it fails to predict the velocity gradient that we observe. It is possible that 
the Globule may be rotating slightly while undergoing an implosion due to a non-symmetric 
increase in external pressure. More detailed modeling with the inclusion of a non-symmetric 
compression wave is necessary to test this idea. 



4-3.2. The merging of two subsonic flows? 

Although we cannot definitively rule out that the ring may have formed via an external 
compression wave, the above model in its current form is unable to completely characterize 
the observed kinematics within this Globule. Because the extinction ring correlates with the 
C^^O velocity dispersion, and, thus, appears organized, it is unlikely that the extinction ring 
is simply a chance super-position of material along the line of sight. 

Rather than the implosion of the Globule due to an increase in its external pressure, it 
appears that we may be seeing the shearing, or merging, of two separate velocity features. 
Gaussian fits to the spectra across the Globule suggests there are two components that arise 
within the ring; one at ~ —5.7 kms~^ and the other at ~ —5.9 kms~^ (see Fig. [3]). The 
emission in these channels, shown in Figure [21 reveals that while both of these components 
peak within the ring, they also show emission external to, and on opposing sides of, the ex- 
tinction ring. Because both of these components have non-thermal subsonic one- dimensional 
velocity dispersions, we speculate that these components correspond to two distinct subsonic 
flows of turbulent gas that are merging in the center of the Globule. 

Of all the measured parameters, we find the spatial match between the one-dimensional 
velocity dispersion and the visual extinction to be the best. Thus, the convergence of these 
two flows within the extinction ring has potentially caused the enhanced column density. It 
may be, therefore, that the extinction ring is neither a ring nor shell seen in projection, but 
simply arises due to the merging of molecular gas and dust as the two flows are interacting. 
Thus, the increase in the velocity dispersion and non-thermal pressure toward the extinction 
ring may simply reflect the fact that two flows are interacting here. Because the outer regions 
of the Globule appear thermally dominated and stable against collapse, it may be that the 
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Globule is evolving from the outside in and has yet to stabilize, as evidenced by the lack of 
a centrally condensed core. 

If the two flows were indeed interacting within the ring, then we would also expect to see 
an increase in the C^^O integrated intensity toward the ring: in the absence of temperature, 
opacity and depletion effects, higher velocity dispersions should translate directly into higher 
integrated intensities. Although a general increase in the integrated intensity is clearly seen 
toward the ring, it is interesting to note that the C^^O emission and extinction do not peak 
at the same location within the ring. Figures [5] and [7] show that the C^^O and CS integrated 
emission peaks toward western arc of the ring whereas the extinction and velocity dispersion 
peak toward the eastern arc. These morphological differences can be explained in terms of 
molecular depletion at high extinctions. Figure [S] shows that both the C^^O and CS are 
depleted for Ay > 9.5 mags. Because the eastern arc is dominated by extinctions > 9.5 
mags it shows a higher degree of molecular depletion compared to the western arc. 

As mentioned previously, there is a temperature gradient across the Globule that goes 
from a Tk of 9.6 ± 0.5 K in the north-west to ~ 7.6 ± 0.4 K in the south-east. Thus, 
it appears that the two distinct velocity flows may have slightly different temperatures. 
Moreover, the fact that the temperature does not increase significantly within the ring, 
supports the idea that it has not formed due to the propagation of a shock front through 
the Globule. If this were the case, we would expect to see the temperature increase in the 
shock front and be coincident with the ring. 



Conclusions 



Globule 2 is one of the densest cores within the Coalsack Nebula, but contains no 
evidence for star-formation. Despite the fact that overall it appears stable against collapse, 
a recent high angular resolution extinction map reveals Globule 2 contains a central ring of 
enhanced column density and, thus, may in fact be a dynamically unstable, transient core 
on the verge of condensing to form a protostar. It has been speculated that this central ring 
of column density has forme d either as a result of turbulence or frorn a shock induced via an 
external compression wave (ILada et al.l |2004J : iHennebelle et al.l |2006| ) . 



Here we have combined dust extinction measurements with a high angular resolution 
q18q molecular line map in order to study the detailed kinematics of Globule 2. The 
q18q j-Qa^p reveals that the Globule is marginally gravitationally bound and not virialized. 
Moreover, the C^^^O map reveals a complex velocity profile that appears to arise from two 
distinct velocity features that are present both toward the ring and external to it. In addition. 
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these velocity features appear to trace material on opposing sides of the extinction ring. 

Each of these velocity flows have non-thermal subsonic velocity dispersions. Outside 
the extinction ring the Globule is thermally dominated, however, both the measured velocity 
dispersion and non-thermal pressure increase toward the extinction ring. Of all the measured 
properties, we flnd the best spatial correspondence between the one-dimensional velocity 
dispersion and the visual extinction. We suggest, therefore, that the extinction ring is an 
enhancement in column density due to the merging of the two subsonic flows. Thus, it 
appears as though we are seeing the convergence of two subsonic flows within a bound 
object. Moreover, because the outer regions of the Globule appear stable against collapse 
and there is no centrally condensed core, it may be that the Globule is evolving from the 
outside in and has yet to stabilize. 

Finally, although we cannot entirely rule out recen t dynamical models for the formation 



of the extinction ring via an external compression wave (IHennebelle et al.ll2006l ). these models 
are unable to completely characterize the observed C^^O kinematics. The discrepancies 
between the model and the data may arise because of the simplifled assumption of a spherical 
compression wave. 

We thank the referee. Dr. Whitworth, for a thorough reading of the paper and the useful 
suggestions which have improved the paper considerably. We also thank Alyssa Goodman 
and Scott Schnee for providing us with the IDL gradient fltting code. 
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Right Ascension (J2000) 

Fig. 1. — The visual extincti on (Ay) map of Globule 2 in the Coalsack Nebula (15" angular 
resolution; iLada et al.l |2004| ) in color scale overlaid with contours of the C^^O integrated 
intensity (46" angular resolution). The emission was integrated over the velocity range of 
—7.0 to —5.0 kms~^. The contour levels are from 1.0 to 2.0 in steps of 0.1 Kkms~^. Note 
that the general morphology of the C^^O integrated intensity follows the enhanced extinction 
associated with the ring. 
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Fig. 2. — C^^O (1-0) channel maps overlaid with contours of the smoothed visual extinction 
map toward Globule 2 (46" angular resolution). The contour levels are from 6 to 12 in steps 
of 0.75 magnitudes. The color maps show a single channel centered at the vlsr marked 
above each image. All images cover the same extent in Right Ascension and Declination (~ 
6' X 6') and have the same color scale from a T^b of (black) to 3.4 K (pink). We see a 
clear velocity gradient across the Globule. 
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Fig. 3. — Normalized number distributions for the peak temperature (T^^), central velocity 
{vlsr), and one-dimensional velocity dispersion (cr) output from the Gaussian fitting pro- 
cedure. In all cases the histograms with horizontal lines represent the data for the 'broad' 
line component, while histograms with diagonal lines represent the 'narrow' line component. 
Note that the two components show a range in their central vlsr, suggesting the kinematics 
of the Globule are more complicated than two simple Gaussian emission profiles. 
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Fig. 4. — C^^O (1-0) position- velocity, {i-V), diagram (color-scale and contours). The emis- 
sion was averaged over the extinction ring in both Declination (left) and Right Ascension 
(right). A velocity gradient is most apparent in the left panel (the Declination-averaged 
position- velocity diagram) which is not surprising given this direction is roughly perpendic- 
ular to the axis of the apparent gradient in the first moment map (see Fig. [5l). 
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Fig. 5. — C^'^O (1-0) moment maps in color scale with contours of the smoothed extinction 
image. Top left Zeroth moment (integrated intensity). Top right: First moment (intensity 
weighted velocity field). Bottom left: Second moment (velocity dispersion). Bottom right: 
Skewness image. The contour levels are from 6 to 12 in steps of 0.75 magnitudes. The 
extinction ring is seen as an enhancement in the C^^O integrated intensity. The first moment 
map reveals a velocity gradient across the Globule. The velocity dispersion measured by the 
second moment map clearly increases toward the extinction ring. The skewness image reveals 
that the profiles change from a negative (blue) asymmetry to a positive (red) asymmetry 
across the Globule. 
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Fig. 6. — Kinetic temperature (T^) across the Globule derived from the optically thick ^^CO 
emission. Overlaid are contours of the smoothed visual extinction image (levels are from 6 
to 12 in steps of 0.75 magnitudes). We see a temperature gradient across the Globule of 9.6 
± 0.5 K in the north-west to ~ 7.6 ± 0.4 K in the south-east. 
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Fig. 7. — CS integrated intensity image overlaid with contours of the smoothed visual ex- 
tinction image. The contour levels are from 6 to 12 in steps of 0.75 magnitudes. Bright 
CS emission is seen coincident with the extinction ring, however, their overall morphologies 
are different. The morphology of the CS integrated intensity emission matches very well the 
C^^O integrated intensity. 
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Fig. 8. — C^^O (left) and CS (right) integrated intensity as a function of visual extinction 
(Ay) for all positions across the Globule. The individual data points are shown as small 
dots while the filled circles were generated by taking the median of the integrated intensities 
within evenly-spaced Ay bins (0.7 mags). The errors bars in Ay represent the range in each 
bin. The error bars in the integrated intensity represent the dispersion in the data. The 
dashed vertical line marks the approximate visual extinction at the edge of the extinction 
ring (~ 9.5 mags). The solid line is a least-squares fit to the data for 4 mags < Ay < 9.5 
mags. We see a general correlation between the C^^O and CS integrated intensity and Ay 
over the complete range in Ay, however, we also note that both the C^^O and CS appear 
depleted at high extinctions (Ay > 9.5 mags). 
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Fig. 9. — Visual extinction (Ay), C^^O integrated intensity, and the ratio of C^'^^O integrated 
intensity to Ky across the Globule. These plots were generated by taking the mean values 
in a 45" strip of constant Declination across the Globule. The Ay across the Globule 
shows the enhanced extinction toward the ring and its central 10% decrement. The C^^O 
integrated intensity does not show a similar decrement. The lower panel shows that the ratio 
is essentially constant toward the western arc of the ring (positive offsets), but drops by ~ 
20% toward the eastern arc (negative offsets). 
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Fig. 10. — Measured non-thermal velocity dispersion {ant', upper panels) and the ratio of 
thermal to non-thermal pressures (Rp; lower panels) as a function of Ay for all positions with 
the Globule (C^^O left, CS right). The individual data points are shown as small dots while 
the filled circles were generated by taking the median of the integrated intensities within 
evenly-spaced Ay bins. Each Ay bin is 0.7 mag wide; the errors bars in Ay represent the 
range in each bin. The error bars in a^t and Rp represent the dispersion in the data. In the 
upper panels the dotted horizontal line marks the one-dimensional isothermal sound speed 
(0.19 kms~^ for hydrogen in a 10 K gas). In the lower panels the dotted horizontal line marks 
the point at which thermal and turbulent pressure are equal. The dashed vertical hne marks 
the approximate visual extinction at the edge of the extinction ring (~ 9.5 mags). For C^^O, 
we see that as Ay increases the non-thermal velocity dispersion increases while the ratio of 
thermal to non-thermal pressure decreases. Coincident with the ring (Jnt transitions from 
subsonic to trans-sonic and Rp transitions from greater than to less than 1. This suggests 
that thermal pressure dominates external to the extinction ring and that turbulence may be 
more important within the ring. In contrast, ant and Rp appear constant across the Globule 
as traced by the CS emission. 
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Fig. 11. — Comparison between our C^'^O data and the models of iHennebelle et al.l (120061 ) 
for different impact parameters, b, across the Globule. The spectra (right panels) were 
generated by averaging over annuli at the impact parameters 6=0, 0.04, 0.08 and 0.12pc 
(shown as circles on the extinction image : left panel). The example model spectra were 
estimated from figure 2 of IHennebelle et al.l (120061 ) and are plotted so that the center velocity 
(their V=0.0 kms~^) corresponds to the v^^^ of the Globule (—5.8 kms~^). The individual 
line centers, widths, and relative intensities were taken from their figure, while the peak 
temperatures were scaled to approximately match the data. We find that at larger impact 
parameters (6 > 0.08 pc) the model matches the data well. Toward the ring (6 < 0.04pc), 
however, we see noticeable differences between the data and the model which might arise from 
the initial assumptions of the simple model (i.e. inflow speeds and a symmetric compression 
wave). 



